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Abstract 

Hyperbolic structures are obtained by tiling a hyperbolic surface with negative Gaussian curva- 
ture. These structures generally exhibit two percolation transitions: a system-wide connection can 
be established at a certain occupation probability p = p c \ and there emerges a unique giant cluster 
at Pc2 > Pel- There have been debates about locating the upper transition point of a prototypical 
hyperbolic structure called the enhanced binary tree (EBT), which is constructed by adding loops 
to a binary tree. This work presents its lower bound as p C 2 > 0.55 by using phenomenological 
renormalization-group methods and discusses some solvable models related to the EBT. 
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I. INTRODUCTION 



Percolation has been one of the most popular model systems in statistical physics and 
it still remains as an active research area. For a classical introduction, one may refer to 
Refs. l|, |2[ . One recent observation in this field is that there generally occur two percolation 
transitions if the size of a given system expands exponentially fast as its length scale grows. 
For example, the size of a binary tree increases as iV ~ 2 L with the number of layers L 
[Fig. CE^a)]. At the first percolation point p^ 66 = 1/2, it becomes possible to establish a 
global connection, and the resulting cluster size scales linearly with L. However, this cluster 
size is still negligible compared with N since L/2 L — > as L increases. We find the largest 
cluster size Si comparable to N only at p = 1, which determines p^ 2 ee = 1. This tree in fact 
belongs to a category called hyperbolic lattices, obtained by tessellating a hyperbolic surface 
appearing in hyperbolic geometry. Since such double percolation transitions in hyperbolic 
structures were revealed by numerical calculations [3] , there have been debates about locating 
the upper transition point (4 - T] , particularly by dealing with a prototypical lattice model 
called the enhanced binary tree (EBT). This structure is not a tree in itself, but is derived 
from the binary tree by connecting vertices on the same layer horizontally [Fig. QJb)]. It thus 
describes spreading along a branching structure with possible horizontal transfer. While the 
duality relation implies p c2 = 0.564(1) [4], we have obtained p C 2 ~ 0.5 by utilizing a simple 
extrapolation of the largest cluster size Si ~ N~^, which is correct for a tree This 
looks also consistent with the observation of s 2 /s\ where S2 is the size of the second largest 
cluster (5)]. We have even tried to explain this estimate p c2 = 1/2 analytically in combination 

7| , pointing 



with numerical observations and approximate renormalization-group methods 
out that the duality argument does not have a solid mathematical ground here. 

Recently, Ref. js| revisited this issue by calculating the crossing probability. According 
to the conformal field theory {9], the crossing probability for a unit disk whose boundary is 
divided at four points Zi,z 2 , z 3 and Z4 is given by 

A — £(§) — i nFif- -■ - n) (i) 

r(|)r(|)^ 2 ll 3'3'3'^' { > 

with the gamma function T, the hypergeometric function 2 F±, and the cross ratio 

(21 - z 2 ){z z - z 4 ) 



n 



(zi - z 3 )(z 2 - Z A ) ' 



2 



FIG. 1: (a) Schematic representations of a simple binary tree and (b) of the EBT derived from (a). 
These are drawn on the Poincare disk, where the circular boundary indicates points at infinity, (c) 
A part of the EBT. 



If the boundary is divided into four equal pieces, e.g., Z\ = —1, z 2 = —i, z 3 = 1, and Zj 



the cross ratio becomes r\ = 1/2 and we immediately find R = 1/2 (see, e.g., Refs. [lOj, 111] 
on the connection between the hyperbolic geometry and the conformal field theory). Such a 
point where R = 1/2 is denoted as a duality point in Ref. 8|. For each of several hyperbolic 
structures considered there, they have numerically calculated R(p) by dividing the boundary 
into four equal intervals. Then by extrapolating R(p) to the large-size limit at the inflection 
point, Ref. [8j] suggests that the limiting tangent line gives an upper bound of p cl and a 
lower bound of p c2 . A notable point is that the slope of the line converges to a finite value, 
which clearly differs from the two-dimensional (2D) results. This method yields p c2 > 0.503 
for the EBT, questioning the validity of the claim that p C 2 = 1/2. They have also estimated 



. — e with e < 1 as growing 
4j. In this paper, equipped 



p C 2 as 0.564(10) by extrapolating the value of p where R(p) = 
the system size, which is consistent with the estimate in Ref. 
with better analytic tools than before, we too reach a conclusion that p c2 is indeed larger 
than 1/2. Our new lower bound, p c2 > 0.55, is obtained by transfer-matrix calculations for 
percolation and includes the lower bound in Ref. [8(. 

This paper is organized as follows. In Sec. Ull we consider two solvable models. Even 
though both of them have the trivial transition point p c2 = 1, this consideration gives an 
insight about percolation in the EBT. Then in Sec. IHIl we deal with the EBT in two different 
ways: one is the block-cell transformation and the other is the transfer-matrix method. Both 
of them lead to p c2 > 1/2 but the latter gives a sharper bound. We then conclude this work 
by reexamining our previous estimate in Sec. HVl 
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FIG. 2: (Color online) (a) Schematic representation of a ternary tree, (b) Crossing probability R 
given by Eq. ([2|) with a tangent line at the inflection point p ~ 0.524 032 (dotted line). 

II. SOLVABLE MODELS 



A. Ternary tree 

A tree with coordination number g = 4 is the simplest example to calculate the crossing 
probability. The coordination number is chosen to provide the structure with natural four- 
fold symmetry. The center node has four branches [solid lines in Fig. EJ^a)], each of which 
leads to a tree with branching ratio b = g — 1 = 3 [dotted lines in Fig. EJ^a)]. Inside each 
tree, the probability ip to connect the top to the boundary can be described by the Galton- 
Watson process 12], where the extinction probability w can be identified with 1 — ip. The 
number of offsprings k for each node is chosen from a binomial distribution -6(3, p), where 
p is the occupation probability of each bond. The generating function is readily obtained as 



(s) = X]fc^=o sk (l)p k (^ ~ p) 3 ^ = (1 — P + sp) 3 . The extinction probability is then given by 



a solution of the equation 4>{w) = w [12], which is 



w 



or, equivalently, 



if < p < | 

_3p2 +2 p3 + v / 4p 3_ 3p 4 1 

2p3 11 3 - P - i; 



if < p < | 

2p3 11 3 - P - L - 



If we also consider the probability to connect two opposite branches of the central cross in 
Fig. E^a), we get the crossing probability as 



R(p) = 2p 2 (l - py^p 2 + Ap A (l - p)ip 2 + p 4 [2^ 2 (l - V) + 4 ^ (1 - tp) + n 



(2) 
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FIG. 3: (Color online) Binary tree with a ring at the boundary. (a),(b) Probabilities of connection 
in A ABC determine those in AA'B'C', resulting in the recursion Eq. ([3]). (c) The probability of 
connection between two farthest points on the boundary, P(BC), when the recursion equation is 
iterated n times, (d) The same quantity for p > 0.9 with n = 4, . . . , 7. 

where 2p 2 (l — p) 2 , 4p 3 (l — p), and p 4 describe the connecting configurations of the central 
cross and the other ^-dependent parts describe configurations of the trees attached to the 
cross. One should note that we have considered crossing in either direction, while it is in 



only one given direction in Eq. ([I]) and Ref. 



This difference in the definition of crossing 



will not change any essential behavior, however. We plot the result in Fig. MJd), and an 
interesting point is that the slope of this function is finite everywhere between p cl = 1/3 and 
p C 2 = 1, in accordance with the numerical analysis in Ref. 8]. Note that this is markedly 
different from the 2D percolation where the slope diverges at the critical point. We also see 
that the tangent line at the inflection point p « 0.524 032 does give a lower bound of p C 2 as 
well as an upper bound of p c \ as suggested in Ref. 



B. Binary tree with a ring at the boundary 

Let us add loops to a tree by attaching a ring along the boundary points. This may be 
regarded as a first step toward making the EBT, and even this single ring can introduce a 
large number of loops into the system. We first focus on the smallest triangle touching the 
boundary and denote it as AABC [Fig. 02(a)]. Following Ref. [13|, we define P(ABC) as 
the probability that all the three points are connected to one another, while P(ABC) as 
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the probability that there is no connection among them. In addition, P(ABC) means the 
probability that B and C are connected but A is not. One can also define P(ABC) and 
P(ABC) in the same way. These cover the whole possibilities by 

P(ABC) + P(ABC) + P(ABC) + P(ABC) + P(ABC) = 1. 

Considering a larger triangle AA'B'C containing the two smallest triangles [Fig. [3(b)], we 
find that it is possible to express the five probabilities of AA'B'C with those of AABC. 
Note that the left-right symmetry is preserved by this transformation so that we have three 
independent variables x = P(ABC), y = P(ABC), and z = P(ABC). After some algebra, 
the transformation turns out to be 

/ 2 / 2 2 2 \ 

x = p (z + 2xz + 2xy — 2px + 3x ), 

y' = p{y-px + xf } (3) 
z = —p(pz 2 — pyz + p 2 xz + pxz — z + pxy — p 2 x 2 + 2px 2 — x), 

where the prime is in order to indicate probabilities for AA'B'C. The initial condition is 
given by counting the possibilities in AABC as 

x = p 3 + 3p 2 (l — p), 
y = p(l-p) 2 , 
z = p(l — p) 2 . 

The quantity of interest is P(BC) = P(ABC) + P(ABC) = x + y, and this can be obtained 
exactly at every iteration step [Fig. W^c)}. Note that this quantity is closely related to the 
crossing probability since it measures the chance for a boundary point to connect to another 
boundary point far away, which is possibly achieved through the inner part of the system. 
When this transformation is iterated, we observe that P(BC) eventually vanishes except at 
p = 1 [Fig. E(c)], so we conclude that the added loops are not enough to make p C 2 nontrivial. 
However, it is notable that the convergence is so slow that it is hard to determine p c i by 
naive extrapolation [Fig. [3(d)]. 
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FIG. 4: (Color online) (a) Block-cell transformation for eight bonds, where r n is the renormalized 
connection probability at the n-th iteration step, represented by double lines, (b) Transformation 
for 23 bonds and (c) its application at p = 0.4 and 0.5. The dotted line indicates a slope of 1. (d) 
A zoomed view at p = 0.5 shows that this p is still below p C 2- 

III. ENHANCED BINARY TREE 
A. Block-cell transformation 

The block-cell transformation shown in Fig. IH(a) was introduced to get a lower bound of 



p C 2 in Ref. 



71 



It yields a lower bound because we systematically overestimate connection 



at each transformation. Using this transformation in Fig. H^a), we concluded p C 2 > 1/2 fl, 
because the limiting connection probability became one for p > 1/2. Later in Ref. [14J, it 
was pointed out that the same method was applicable to find lower bounds for the usual 2D 
percolation thresholds as well, and also that such a lower bound approached the true critical 
point p c as we used a larger block. For the square lattice, a small block already predicts the 
correct critical point of the bond percolation, p^ uarc = 1/2, and using a larger block does 
not change this estimate, which can be an indication of the exactness of p^ uarc = 1/2. This 
observation motivates us to take a larger block in the EBT [Fig. H](b)], which requires us 
to check 2 23 « 8 x 10 6 configurations. The enumeration is straightforward with a personal 
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computer and the result is as follows: 



r n+1 = 9p l2 r 7 n - 91p u r 7 n + A09p l0 r 7 n - W71p 9 r 7 n + 1795p s r 7 n - 1982p 7 r 7 n + UUp 6 r 7 n 
-590p 5 r 7 n + 9\p 4 r 7 n + 25p 3 r 7 n - 7p 2 r 7 n - 3pr 7 n + r 7 n - 25p 12 r^ + 227p u r 6 n 
-90V°r£ + 2060p 9 r£ - 2928p 8 r£ + 2632p 7 r 6 n - UUp 6 r 6 n + 35Ap 5 r 6 n + 25p 4 r 6 n 
-27p 3 r 6 n - 2p 2 r 6 n + 2pr 6 n + 2Ap 12 r 5 n - 185p n r 5 n + 600p 10 r^ - 1038p 9 r 5 n + 969p 8 r 5 n 
-35Ap 7 r 5 n - 165p 6 r* + 210p 5 r^ - 57p 4 r 5 n - 9p 3 r 5 n + 5p 2 r 5 n - 9p 12 r 4 + 5lp ll r 4 
-102p 10 r 4 + 60p 9 r 4 + 60 An - 70p 7 r 4 - A6p 6 r 4 + 92p 5 r 4 - 35p 4 r 4 - 5p 3 r 4 
+Ap 2 r 4 n + p 12 r 3 n - 2p n r 3 n - 2Ap 9 r 3 n + ll5p s r 3 n - 198p 7 r 3 n + 150p 6 r„ - 32p 5 r 3 n 
-I8p 4 r 3 + 8p 3 r 3 - 3p 10 r 2 n + 13p 9 r 2 n - Up % r 2 n - I2p 7 r 2 n + 3lp & r 2 n - I3p 5 r 2 n 
— 6p 4 r^ + Ap 3 r 2 + 3p 8 r n — lQp 7 r n + 31p 6 r n — 2Qp 5 r n + 8p 4 r n — p 6 + 5p 5 
-8p 4 + Ap 3 + p. 



The system-wide connection probability grows as we increase p [Fig. EKc)]. By checking 
the value of p where becomes 1, we locate a lower bound of p c2 . In fact, a careful look 
shows that p = 1/2 is still below p c2 [Fig. H](d)] and locates a sharper bound p c2 > 0.523. 

B. Transfer-matrix method 

In studying the Ising model on the EBT in Ref. [15] ]. we pointed out that the transfer- 
matrix method performed better than the block-cell transformation. Hence, we employ 



the transfer- matrix formalism developed for percolation in Ref. [161 ] . First, we consider a 
unit cell of three spins as shown in Fig. [^a). Note that a bond on the bottom line has 
a different probability of q n from the others with p. By attaching these cells from left to 
right, we construct an indefinitely long strip, or a layer of width 1, which we can solve by 
using the transfer-matrix method. When we consider connection to the leftmost side, this 
cell has three possibilities: first, only the top point A is connected (case 1), second, only 
B is connected (case 2), or finally, both of them are connected (case 3). So we have nine 
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FIG. 5: (Color online) (a) Unit cell for a transfer matrix calculation of a layer with width 1. Each 
bond at the bottom line has an occupation probability q n (represented by a double line) , while the 
other bonds have an occupation probability p. (b) Iteration maps for p = 0.4 and 0.5. (c) p = 0.5 
has an intersection at q^ < 1 and therefore lies below p C 2. The dotted lines indicate q n +\ = q n - 

possibilities of connection in total as follows: 



Pi 


= P(1 


->• 


1) 


= P(p 2 Q 2 a-Pql-PQn + 1), 


P21 


= P(1 


— >• 


2) 


= pg n (l -p)(l + q n -pq n ), 


P31 


= P(1 


— >■ 


3) 


= p 2 q n (l + q n -pq n ), 


P12 


EE P(2 


— >■ 


1) 


= p 2 (l- g„)(l + 2q n - 2pq n ), 


P22 


EE P(2 


-)• 


2) 


= QnO- - p){p 2 + q n + pq n - 2p 2 q n ), 


P32 


EE P(2 


->• 


3) 


= P 2 q n {p + 2g„ - 2pq n ), 


Pl3 


EE P(3 


— >■ 


1) 


= P0--Qn)0- + Qn-PQn), 


P23 


EE P(3 


->• 


2) 


= qn(i-p)(p + q n -pqn), 


P33 


EE P(3 


->• 


3) 


= pq n {p + qn-pq n )- 



The global probability of connection to the leftmost side when n(^> 1) blocks are attached 
will behave as ~ A™, where A is the largest eigenvalue of this 3x3 matrix {Py}. So we 
replace this layer by a one-dimensional chain, and identify its occupation probability q n +i 
with A to recover the original configuration in Fig. 13(a), but with q n+ i instead of q n . This 
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FIG. 6: (Color online) (a) Unit cell for a transfer-matrix calculation of a layer with width 2. (b) 
Ten possible cases of connection. A black dot means that it is connected to the leftmost side, and 
a white dot means it is not. A double circle means that these two are connected to each other, 
while none of them are connected to the leftmost side, (c) The resulting iteration map suggests a 
new lower bound of p C 2 as around p = 0.55. 

iteration therefore determines q n +i as a function of q n and p, and we will find a limiting 
value goo = Hindoo q n for a large system. This renormalized connection probability is an 
increasing function of p [Fig. GD[b)]. Again, we are interested in the value of p making = 1, 
and such p is found to be ~ 0.504. In short, this confirms that p = 1/2 is strictly below p c2 
[Fig. He)]. 

We can take a larger layer, expecting a sharper bound [Fig. HJa)]. This consideration has 
ten possible listed in Fig. [6](b). Here, the black dots are connected to the leftmost 

side while the white dots are not. It is important to consider possibilities that two white dots 
may be connected to each other since a percolating path may go backward for a while, so such 
dots are represented by double circles. In fact, the case indexed as 4 is not accessible from 
any other states so it can be discarded. Each matrix element is expressed as a high-order 
polynomial, so the average number of terms per polynomial amounts to 18.4. It is merely 
a mechanical procedure to obtain the matrix elements so we show only the final result with 
the largest eigenvalue, which is identified with q n +i as above. The result shown in Fig. [6]^c) 
shows the highest lower bound of p C 2, which is around p = 0.55. It includes the numerical 
lower bound suggested in Ref. 8] and the estimate based on the duality relation J4]. 
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FIG. 7: (Color online) Ratio of the second largest cluster size with respect to the largest cluster 
size in the EBT with different numbers of layers L. Our calculation suggests that the slope at 
p = 1/2 will converge to a finite value in the large- ./V limit. 



IV. DISCUSSION 



- S3, 



At the time of writing Refs. |3|, |5[, we assumed that the percolating properties in hyper- 
bolic lattices could be inferred from known 2D results and also from the results of a tree. 
For example, S2/S1 has a diverging slope at the emergence of a giant cluster both for a 2D 
plane and for a tree. That is why we expected the same behavior for hyperbolic structures 
as well. However, as we see a clear difference from the 2D result in Fig. E^b), such an 
assumption now looks quite dubious. Based on the results so far obtained, it seems more 
plausible that this ratio also has a constant slope at the large- N limit (Fig. [7]). This implies 
that percolation in the EBT is neither similar to its 2D counterpart nor to the percolation in 
a simple tree. In particular, we see that the competition between the largest and the second 
largest clusters appears milder than has been believed, so that S2/S1 may vanish smoothly 
around p C 2- In addition, it is inevitable to reconsider the phenomenological description of the 
critical phenomena around p C 2 by using scaling collapse jsj since it is likely that the critical 
points in the hyperbolic lattices have been generally underestimated. Our second example 
in Sec. HI1 however, suggests that it can be difficult to extract the critical behavior if one 
solely relies on numerical data. It will be interesting to challenge this problem by making 



use of the recent analytic approaches to hierarchical structures (see, e.g., Ref. [171. Il8|). 
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